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ABSTRACT 

A method to compute several scalar quantities of Cosmic Microwave Background 
maps on the sphere is presented. We consider here four type of scalars: the Hessian 
matrix scalars, the distortion scalars, the gradient related scalars and the curvature 
scalars. Such quantities are obtained directly from the spherical harmonic coefficients 
aim of the map. We also study the probability density function of these quantities 
for the case of a homogeneous and isotropic Gaussian field, which are functions of 
the power spectrum of the initial field. From these scalars it is posible to construct 
a new set of scalars which are independent of the power spectrum of the field. We 
test our results using simulations and find a good agreement between the theoretical 
probability density functions and those obtained from simulations. Therefore, these 
quantities are proposed to investigate the presence of non-Gaussian features in CMB 
maps. Finally, we show how to compute the scalars in presence of anisotropic noise 
and realistic masks. 



Key words: 

ground 



methods; analytical - methods: statistical - cosmic microwave back- 



1 INTRODUCTION 

The Cosmic Microwave Background (CMB) is currently one 
of the most valuable tools of cosmology, providing with a 
great wealth of information about the universe. A particu- 
larly interesting subject is whether the CMB temperature 
fluctuations follow or not a Gaussian distribution. This is 
a key issue since Gaussian fluctuations are predicted by 
the standard inflationary theory, whereas alternative the- 
ories produce non-Gaussian signatures in the CMB. In ad- 
dition, foregrounds and systematics may also introduce non- 
Gaussianity, which should be carefully studied in order to 
avoid its misidentification with intrinsic non-Gaussian fluc- 
tuations. 

A large effort has been recently devoted to the study 
of the Gaussian character of the CMB using the multi- 
frequency all-sky CMB data provided by the WMAP satel- 
Ute of NASA (Bennett et al. 2003a), flnding, in some 
cases, unexpected results. Some authors have found that the 
WMAP data are consistent with Gaussianity using different 
types of analysis (Komatsu et al. 2003, CoUey & Gott 2003, 
Gaztahaga & Wagg 2003, Gaztariaga et al. 2003). How- 
ever, other studies have produced a positive detection of 
non- Gaussianity and/or have shown north-south asymme- 
tries (Chiang et al. 2003, Eriksen et al. 2004a,b, Park 2004, 
Copi et al. 2004, Vielva et al. 2004, Hansen et al. 2004, 
Mukherjee & Wang 2004, Land & Magueijo 2004, Hansen 



et al. 2004, Larson & Wandeff 2004, McEwen et al. 2004, 
Cruz et al. 2005). Although some of these results could be 
explained by the presence of foreground contamination, in 
other cases the origin of the detection is uncertain and a 
primordial origin can not be discarded. These results mo- 
tivates, even more, the development of novel techniques to 
perform further Gaussianity analyses of the CMB. 

In this paper we have focused on the study of sev- 
eral scalar quantities constructed from the derivatives of the 
CMB field on the sphere. In particular, we have considered 
the modulus of the gradient, the Laplacian, the distortion, 
the shear, the ellipticity, the shape index, the eigenvalues 
of the negative Hessian matrix, the Gaussian curvature and 
the derivative of the squared modulus of the gradient. A 
procedure to calculate the scalars directly from the atm har- 
monic coefficients is provided. In addition, we have obtained 
expressions for the probability distributions of the previous 
quantities for the case of a homogeneous and isotropic Gaus- 
sian random field (HIGRF), which depend on the power 
spectrum of the initial field. In order to remove this de- 
pendence, we have also constructed a new set of normalized 
scalars whose probability distribution is independent of the 
power spectrum. 

A number of works have been already devoted to the 
study of some of these or other related quantities. Also, some 
of them have been applied to study the Gaussianity of the 
WMAP data, providing interesting results. In a pioneering 
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work, Bond & Efstathiou (1987) studied different statistical 
properties ol the CMB, assuming a 2-dimensional Gaussian 
field, which included the number densities of hot and cold 
spots, the eccentricities of peaks and peak correlation prop- 
erties. Barrciro ct al. (1997) studied also the mean number 
of maxima and the probability distribution of the Gaussian 
curvature and the eccentricity of peaks of the CMB for dif- 
ferent shapes of the power spectra. Wandelt, Hivon & Gorski 
(1998) developed efficient algorithms for fast extrema search 
as well as for the simulation of the gradient vector and the 
curvature tensor fields on the sphere associated to the tem- 
perature field. The power to detect non-Gaussianity in the 
CMB of the number, eccentricity and Gaussian curvature 
of excursion sets above (and below) a threshold was tested 
by Barreiro, Martmez-Gonzalez & Sanz (2001) using Gaus- 
sian and non-Gaussian simulations, finding that the Gaus- 
sian curvature was the best discriminator. Dore, Colombi 
& Bouchct (2003) tested the power of a technique based 
on the proportion of hill, lake and saddle points (which are 
defined attending to their local curvature) on flat patches 
of the sky. A study of the ellipticity of the CMB was per- 
formed for the COBE-DMR (Gurzadyan & Torres 1997) and 
Boomerang data (Gurzadyan ct al. 2003), which seemed to 
indicate a slight excess of ellipticity for the largest spots 
with respect to what was expected in the standard model. 
These results were recently confirmed by repeating the anal- 
ysis on the WMAP data using the same region of the sky 
observed by Boomerang (Gurzadyan et al. 2004). Eriksen 
et al. (2004b) applied the Minkowski functionals and the 
length of the skeleton (Novikov, Colombi & Dore 2003) to 
the WMAP data finding evidence of non-Gaussianity and 
asymmetry between the northern and southern hemispheres. 
The length of the skeleton is a quantity defined in terms of 
the first and second derivatives of the field and the previ- 
ous work provided with an algorithm to calculate it on the 
sphere for a CMB map in HEALPix pixclization. Hansen et 
al. (2004) found north/south asymmetries on the WMAP 
data using the local curvature, which were consistent with 
the results of Eriksen et al. (2004b). Finally, Cabella et al. 
(2004) used a method based on the local curvature to con- 
strain the value of the non-linear coupling constant /nl from 
the WMAP data. 

The paper is organised as follows. In §2 we introduce 
the different considered scalars and describe how to calcu- 
late them from the covariant derivatives of a 2D field and, 
in particular, for spherical coordinates. In §3 we obtain an- 
alytical (or semi-analytical) expressions for the probability 
density function (pdf) of the considered scalars for a homo- 
geneous and isotropic Gaussian field on the sphere. In §4 we 
construct a new set of related scalars which arc independent 
of the power spectrum of the field and we obtain their theo- 
retical distribution functions for the Gaussian case. The the- 
oretical results are also compared with CMB simulations. §5 
shows how to deal with more realistic simulations, includ- 
ing anisotropic noise and a mask that covers the Galactic 
plane and the point sources. In §6 we present our conclusions 
and outline future applications of this work. Appendix A 
gives the expressions necessary to calculate the derivatives 
of the field from the aem coefficients. In Appendix B we 
present some useful results regarding spherical harmonic se- 
ries. Finally, Appendix C gives some guidelines on how to 



deduce the distribution functions of ordinary and normalised 
scalars. 



2 SCALARS ON THE SPHERE 

Let us consider a 2-dimensional field T(a;i,X2). From the 
field derivatives we cam construct quantities that are scalars 
under a change of the coordinate system (i. e. regular general 
transformation (xi,X2) (x'i,x'2): s'(x'i,x'2) = s{xi,X2))- 
Regarding 1st derivatives, a single scalar can be constructed 
in terms of the ordinary derivative Ti. Regarding 2nd co- 
variant derivatives, T.ij, of T(.ti,.T2) can be expressed as a 
function of the ordinary derivatives Tij and the Christoffel 
symbols P^^- as follows 



- ij-' 



(1) 



To construct linear scalars we need to contract the indices 
of these covariant tensorial quantities. 

The scalars that depend on second derivatives, although 

can be expressed as functions of the covariant derivatives, 
are usually defined in terms of the eigenvalues Ai, A2, of the 
negative Hessian matrix A of the field T(a;i,a;2) 



i-T-.ij) . 



(2) 



That is, Al and A2 are the (negative) second derivatives 
along the two principal directions. In the following we will 
assume Ai > A2. 

Taking into account the values of Ai, A2, we can distin- 
guish among throe type of points (e. g. Dore et al. 2003): hill 
(both eigenvalues are positive), lake (both are negative) and 
saddle (Ai > 0, A2 < 0). We will also distinguish between 
saddle points with |Ai| > IA2I, that we will call saddle-, 
and those with |Ai| < IA2I, saddlo+ (the sign -|-,- refers to 
the sign of the Laplacian, see below). 

Hereinafter, wo will consider spherical coordinates 
{d, <j)), for which the metric is given by 



ds^ = d9^ + SVC? ( 



which gives the following non-zero Christoffel symbols: 



p</> p</> 



cos 9 



(3) 



(4) 



It is convenient to define the following quantities for the 

spherical case: 



1 dT 
sin 9 d(t> 

_ 

1 d'^T cos 9 dT 



1 d^T 



+ 



cos 6* dT 



(5) 
(6) 

(7) 
(8) 
(9) 



Bm9d9d(p sin^ 6) 90 ' 

We will express all the scalars as a function of these five 
quantities. All of them, and therefore the scalars, can be 
easily computed from the aern. coefficients of the field (see 
Appendix A). 
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In the following, we present some of these scalars as 
function of the covariant derivatives. We will also relate 
them to the former defined quantities for the spherical case. 

Note that, although we are considering the spherical 
case, the expressions given for the scalars as a function of the 
covariant derivatives are valid for any 2-dimensional surface. 
In particular, if observing a small portion of the sky, we can 
assume that we have a flat patch. In this case, the scalars can 
be easily obtained taking into account that all the Christoffel 
symbols are zero. 

2.1 The Hessian matrix scalars 

We include in this section the algebraical quantities related 
to the Hessian matrix: eigenvalues, trace and determinant. 



As a function of the covariant derivatives, the determi- 
nant can be expressed as 




Finally, for the spherical coordinate system, we can rewrite 
it as 



d = rs-f. (18) 
2.2 The distortion scalars 

We consider here the distortion, the shear, the ellipticity 
and the shape index. They are all related with powers of 

the difference of the eigenvalues Ai — A2, so they give us 
information about the distortion of the field. 



2.1.1 The eigenvalues 

The eigenvalues of the negative Hessian matrix are scalars 
of the field T, and they can be expressed as follows 



2Ai = - {r\) + ^{T\Y - 2 (T^VT^^, - ^^^,^^^) , (lo) 
2A2 = - {T'\) - ^ - 2 {T'\T'^. - r\r',) . (ii) 

For the spherical coordinate system, we can rewrite them in 
the following way: 



(12) 
(13) 



2Ai = -{r + s) + ^{r - s)^ + (2t)2 , 

2A2 = - (r + s) - ^J(r- sf + {2tf . 

Note that Ai > A2. Therefore for values of A2 > we have 
hill points whereas values of Ai < correspond to lake 
points. These eigenvalues coincide with the two principal 
curvatures of the surface at extrema points where the first 
derivatives are zero. 



2.1.2 The Laplacian 

The Laplacian A+ is defined as the trace of the Hessian 
matrix. Therefore it can be expressed as a function of the 
eigenvalues: 



-Ai 



(14) 



Note that negative values of A+ correspond to hill or saddle- 
points of the field, whereas lake and saddle+ points corre- 
spond to positive values of the Laplacian. 

The Laplacian on the sphere can also be written as a 
function of the field covariant derivatives or r and s: 



\+=Tr=r + s 



(15) 



2.1.3 The determinant of A 

Another scalar that can be constructed is the determinant 
of the negative Hessian matrix A, which is given by 

d=detA = A1A2 . (16) 

Positive values of d correspond to hill or lake points of the 
field whereas saddle points are given by negative values of 
d. 



2.2.1 The shear 

An interesting scalar related to the eigenvalues is the shear, 
which is defined as 



= i(Ai-A2)^ 



(19) 



As a function of the covariant derivatives we can express it 

as follows 



y 



— TT'* 1 ^ — i TT'' T'-' — T'-' T'* 1 

A I i] 9 L » 3 i jj 



and for the spherical coordinate system, we can rewrite it as 



y=-[r-s) +t 



2.2.2 The distortion 



(21) 



The distortion is defined as the difference of the eigenvalues 
of the negative Hessian matrix: 



A- 



A2 



(22) 



therefore, by construction < A_ < 00. We can express the 
distortion as a function of the covariant derivatives of the 
original field in the following manner 



A- = VF^J - 2 [T-^^ - r^^rVJ (23) 
using the spherical coordinate system, we can rewrite it as 

A_ = ^{r-sf + {2tf. (24) 

2.2.3 The ellipticity 
The ellipticity is defined as 
Ai — A2 



^"2(Ai+A2)- ^^^^ 
It can take values in the whole real domain. It is straight- 
forward to show that e € (~|iO) ^^'^ ^ € ('^i I) ^'^^ ^^'^^ 
and hill points respectively. Regarding to saddle points, the 
ellipticity takes values in the range e £ (51 00) for saddle- 

and e € (—00, — for saddle+. 

As a function of the covariant derivatives of the field we 
can express the ellipticity in the following way: 
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and using r, s and t in spherical coordinates the eUipticity 

is rewritten as follows: 



2{r + s) 



(27) 



2.2.4 The shape Index 

There are also other scalars that can be constructed from 
the ellipticity. One of the most interesting quantities is the 
shape index (Koenderink 1990), defined as 



t = — arctan 

TT 



(28) 



By definition, the shape index is bound t G (—1,1). Hill 
points correspond to values oftG > ^^^^ points to 

(. £ ( I ' 1) ' saddle- points to t G ( ~ 5 > O) ^^'^ finally saddles- 
points are in the range t € (O, |). 



2.3 The gradient related scalars 

We included here the squared modulus of the gradient field, 
and also another quantity related with its derivative. 



2.3.1 The square of the gradient modulus 

The square of the modulus of the gradient is defined as fol- 
lows: 



2.4.1 The Gaussian curvature 

The Gaussian curvature, also called intrinsic curvature, is 
defined as the product of the two principal curvatures of the 
surface defined by the temperature field. It can be expressed 
as a function of the field derivatives: 



KG = 



2 [1 + {T.Xi)f 



(33) 



wo can rewrite it in the spherical coordinate system as fol- 
lows: 



Kg = 



r s-r 



(34) 



where kg G (— cx),cx)) by construction. Note that for cx- 
trcma, the Gaussian curvature coincides with the determi- 
nant of A. This scalar depends on intrinsic properties of the 
surface. 



2.4.2 The extrinsic curvature 

This scalar is defined as the average of the two principal 
curvatures of the field and it can be rewritten as a function 
of previous scalars in the following way: 



2 ^/TTg 



A4 



+ 



(35) 



This scalar gives an idea of how the surface generated by 
the temperature field is embedded in 5R^. Note that in the 
extrema this quantity coincides with the Laplacian except 
by a constant factor. 



VT 



T'^T.i 



(29) 



It gives information about the smoothness of the field. Tak- 
ing into account equations (5) and (6), g can be expressed 
for the spherical coordinate system as 



2 , 2 
g=p +q 



(30) 



2.3.2 Derivative of the square of gradient modulus 

We also consider another non linear scalar, called Dg, which 
is proportional to the derivative of the square of the gradient 
modulus with respect to the arc associated to the integral 
curves of the gradient. As a function of the covariant deriva- 
tives we can write it in the following way: 



Da = 



d 



t5 



(31) 



ds \2- 

In the spherical coordinate system we can express it as 

Dg=rp'' + 2tpq + sq'' , (32) 

where Dg € (— cx),cx)). This quantity is related to the ex- 
trinsic curvature (see below). 

2.4 The curvature Scalars 

We include in this section the intrinsic and extrinsic curva- 
tures. They give us information about the geometry of the 
surface of the temperature field. 



3 HOMOGENEOUS AND ISOTROPIC 
GAUSSIAN RANDOM FIELD ON THE 
SPHERE 

In this section, we will derive the pdf of the scalars defined in 

the preceding section for a HIGRF on the sphere. Under this 
assumption, it can be shown that the quantities {p, q, r, s, t} 
defined in equations (5) to (9), follow also a homogeneous 
and isotropic Gaussian distribution. In order to calculate the 
dispersions and covariances of these quantities, we define the 
moments at as 

2£+l 



47r 



(36) 



where Ct is the power spectrum of the temperature field. 
Making use of the results given in Appendix B, it can be 
shown that: 



< p p > = < q q> 



1 2 



3 2 1 2 
<r r > = < s s > = -(72 — Tfi , 

o 4 



<tt> 



1 2 1 2 



(37) 
(38) 
(39) 



If the initial field has zero mean, {p, q, r, s, t} have also zero 
mean, and their covariances are zero except for the term 
< r s > which is given by: 



1 2 1 2 
<r s >= -f72 -I- -r(Ti = p <r r > 

o 4 



(40) 
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In order to obtain the pdf's of the scalars, it is useful to 
define a new set of Gaussian variables which are uncorrelated 
among them. In particular, we construct R = r + s, S = r — s 
and T = 2t. It is straightforward to show that crfj = and 

j_T . 2 2 12 2 

that (Tg = (Ty = 2C2 ~ '^1 ■ 

As an illustration, in Fig. 1 we show the theoretical 
distribution functions of the variables q and r as well as 
the ones obtained averaging the normalised histograms of 
20 all-sky CMB Gaussian simulations. The simulations have 
been generated using the HEALPix package (Gorski et al. 
1999) for A^side ~ 256 and using the power spectrum given 
by the best-fit model found by the WMAP team (Spergel 
et al. 2003). The C^'s for this model were generated using 
CMBFast (Seljak & Zaldarriaga 1996). The simulated maps 
were smoothed with a Gaussian beam of full width half max- 
imum (FWHM) equal to 2.4 times the pixel size. Note the 
good agreement between the theoretical distribution and the 
numerical results. 

Since we have expressed all the scalars as a function of 
{p, q, r, s, t} and since they follow a simple Gaussian distri- 
bution, we can almost straightforwardly derive the proba- 
bility distribution function of the considered scalars. 

3.1 The Hessian matrix scalars 

3.1.1 The eigenvalues 

Taking into account equation (12), we can obtain the pdf of 
Ai: 

p{\i) = fe -f-fcAie-"^i [1-f erf (ftAi)] , (41) 
where the constants c, /, k and h are given by 

2 ^ _a2_ 

' " ol + crl ' ^ - 72^'' 

fc = ^ci , ft = , (42) 

V2 0-2 

and erf is the error function. For this distribution function 
we obtain < Ai >= ^/Jcys. 

Analogously, we can obtain the pdf for A2 : 

p (A2) = /e "2 - kMe'"^^ [1 - erf {hM)] . (43) 

In this case, we find < A2 ^^/^'^s- Note that there 
is a symmetry between the pdf's of Ai and — A2, such that 
PAi(a:) =px^[-x). 

The map of the eigenvalues Ai and A2, corresponding 
to the CMB Gaussian simulation in Fig. 2 (smoothed with 
a Gaussian beam of FWHM of 8 deg) can be seen in the top 
left and right panel of Fig. 3 respectively. We have chosen 
a low resolution CMB map for a better visualization of the 
different type of points, that with this smoothing appears 
connected, since otherwise the maps of the scalars are dom- 
inated by the small scale and the large scale structure can 
not be appreciated. Note that negative values of Ai, corre- 
sponding to lake points in the original field, form compact 
regions while the rest of pixels, corresponding to other type 
of points, form a filamentary structure similar to a web. In 
a similar way positive values of A2, corresponding to hill 
points, form compact regions surrounded by the filamentary 
structure. 




■0.02 0.02 

q 




-5 



Figure 1. Top panel: the solid line corresponds to the theoretical 
distribution expected for q, the crosses have been obtained aver- 
aging over 20 all-sky CMB simulations and the error bars give the 
dispersion for the same 20 simulations. Bottom panel: the same 
results are shown for the variable r. 




-^.UH.. ^^^^^^^^^^^^^^^^^ -".ll^. 



Figure 2. A Gaussian CMB simulation (in units of K) smoothed 
with a Gaussian beam of FWHM=8° is shown. 
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3.1.2 The Laplacian 

The Laplacian is the addition of two Gaussian variables with 
the same dispersion and therefore it also follows a Gaussian 
distribution with dispersion <T2. 



1 



27rcr2 



exp 



(44) 



The Laplacian is defined in the interval A+ € (—00,00). 
As an example, wc show the Laplacian of a CMB Gaussian 
simulation in Fig. ,3 (bottom left panel). Note that those 
regions with low values of A+ correspond to regions with 
high positive curvature in the initial temperature field, and 
analogously regions with high values of A+ correspond to 
regions with high negative curvature in T. 



3.1.3 The determinant of A 

Using equation (18), we can construct a semi-analytical ex- 
pression for the pdf of the determinant: 



p (d) = Cde" 



\x'^ -\- d\\ dx 



'0 \P 
where d £ (—00, 00) and the constants are 
2 1 



(45) 



7ri/27r (1 - p2) a?at 
P 1 



1 



a = w — 



1 

2^ ■ 



(46) 



and Ko is the zero order modified Bessel function of second 
kind. 

A map of d for the Gaussian CMB sinnilation shown 
in Fig. 2, is given in Fig. 3 (bottom right panel). Positive 
values of d correspond to hill and lake points, so the most 
curved hills and valleys in the original field are given by the 
highest values of d. 

3.2 The distortion scalars 

3.2.1 The shear 

The shear is the addition of two independent squared Gaus- 
sian variables with the same dispersion. Therefore, it follows 
a xi distribution function with mean and dispersion equal 

to < y >= I 



1 ^2 
4(^2 



1^2 



Piv) = — e" 



(47) 



where y £ [0, 00). The map of the shear of a Gaussian CMB 
simulation is shown in the top left panel of Fig. 4. In the con- 
sidered example, we find that large values of y are concen- 
trated in regions which usually correspond to saddle points 
with high values of | A, | in the original map. 



where A_ is positive by construction. In the top right panel 
of Fig. 4 we show the distortion map for the Gaussian CMB 
simulation of Fig. 2. The physical information contained in 
the distortion is basically the same as that of the shear. In 
fact, both maps show quite similar structure. 

3.2.3 The ellipticity 

Taking into account equation (27), we can obtain the pdf of 
the ellipticity: 



p(e) = 2CTsaij|e| 



i 



+ 4a2 e2 



(49) 



where e G (—00, 00). The map of the ellipticity computed for 
a Gaussian CMB simulation is given in Fig. 4 (bottom left 
panel). It is seen that the largest values of | e | correspond 
to saddle points surrounding hill or lake points which form 
compact regions. 

3.2.4 The shape index 

The shape index pdf can be easily derived from the one of 

the ellipticity: 

cos^d^)^"" 



P(0 





cos (ft) 


4 


sin^(ft)| 



2 , 2 



sin 



(50) 



Note that the shape index is bound and takes values in 
the range t G [—1, 1]. The structure of the shape index for 
a Gaussian CMB simulation is shown in the bottom right 
panel of Fig. 4. This scalar presents a similar structure as the 
one of the ellipticity. Higher values of t correspond to lake 
points in the original temperature field while lower values of 
i correspond to hill points. 

3.3 The gradient related scalars 

3. 3. 1 The square of the gradient modulus 

Taking into account equation (29) , we see that g is given by 

the addition of two independent squared Gaussian variables 

with the same dispersion. Therefore this scalar follows a xi 

distribution with mean and dispersion equal to < g >= 
2 



(51) 



where g G [0, 00). As an illustration of the structure of the 
squared modulus of the gradient, we show the gradient map 
of the low resolution CMB Gaussian simulation of Fig. 2 in 
the top panel of Fig. 5. Note that those regions where the 
temperature field changes rapidly correspond to high values 
of g. 



3.2.2 The distortion 



3.3.2 Derivative of the square of gradient modulus 



The distortion is proportional to the square root of the shear, 
thus its probability density function can be easily obtained 
from the one of the shear. 



piX_) = ^^e- 



(48) 



Dg involves the five Gaussian variables (5) to (9) and its 
probability distribution can be expressed as a function of 
one integral: 



piDg) = CD. 



Jo 



dz 



-exp 



2al 



+ 



(52) 
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Figure 3. Hessian matrix scalars associated to tlie temperature map of Fig. 2 are sliown: the eigenvalues Ai (top left) and A2 (top right), 
the laplacian (bottom left) and the determinant (bottom right). 



'Oo, cxj) and the constant Cug is given by 



2-K O-rCr^ 



(53) 



The structure of this scalar is shown in the map of Fig. 5 
(bottom panel) obtained for a Gaussian CMB simulation. 
Note that this scalar was defined as the derivative of the 
squared modulus of the gradient vectorial field, except by 
a constant factor. In this way it gives us extra information 
about the directional variation of the gradient, which is not 
given by g. 



3.4-1 Gaussian curvature 

The pdf of the Gaussian curvature can not be obtained in 
an analytical form but it can be written as a function of two 
integrals: 



(54) 



where kq G (—00, 00), and I (a) is the following integral: 



I (a) 



dxe 



''Ko{\x^ + a\) 



(55) 



3.4 The curvature scalars 

As shown in equations (34) and (35) the curvature scalars 
are a function of 1+ | VT p, which depend on the units of 
the temperature field and, therefore the shape of the curva- 
ture scalars pdf 's will also non-trivially depend on the these 
units (note that for the previously studied scalars, the units 
of the field enter just as a normalization factor in the pdf). 
For this reason we use adimensional fields to calculate the 
curvatures. In particular, we choose to normalise the ini- 
tial temperature field to unit dispersion to enhance possible 
deviations from Gaussianity. 



The constants are given by 



h = 



1 

2^ " 



P 1 



2 2 

(Tr O't CTp 



1 



(56) 



(57) 



We present in the top panel of Fig. 6 the Gaussian curvature 
corresponding to the CMB simulation of Fig. 2 (which has 
been normalised to unit dispersion, as mentioned before). 
Note that regions with positive values of kq correspond to 
hill or lake points while negative values correspond to saddle 
points. This behaviour is similar to the one found for the de- 
terminant, but for Kg the gradient weights this information 
in those regions where the later takes large values. 
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Figure 4. The distortion scalars associated to the temperature map presented in Fig. 2 are plotted. They are the shear (top left), the 
distortion (top right), the eUipticity (bottom left) and the shape index (bottom right). 



3.4-2 The extrinsic curvature 

The extrinsic curvature is defined by equation (35). To ob- 
tain an analytical expression for the pdf of this scalar is very 
complicated and, therefore, we do not include it here. 

However, we show in the bottom panel of Fig. 6 the ex- 
trinsic curvature map corresponding to the normalised CMB 
simulation of Fig. 2. 



4 NORMALISED SCALARS 

In order to enhance possible non-Gaussian signatures in the 
temperature field, it would be desirable to study quantities 
which are independent of the power spectrum of the anal- 
ysed field. With this aim, we have constructed a new set 
of scalars that have this property and that are related to 
the physical scalars of the previous sections. In addition, as 
will be shown later, these new quantities are mathematically 
simpler, and allows one to deal in a straightforward manner 
with anisotropic fields. 

In some cases the relation between the normalised and 
the ordinary scalar, is only a constant factor. In this case the 
distribution function of these new quantities, can be deduced 
straightforwardly from the old ones. We can include in this 
group the square modulus of the gradient and its derivative, 
the Laplacian, the distortion, the shear and the eUipticity. 
The rest of the new scalars should be defined using other 



normalised scalars, and their pdf 's constructed accordingly 
(see Apendix C). 

In the next subsection, we introduce the normalised 
scalars and give their distribution function for the Gaussian 
case. We have summarized this information in table 1. 

In order to test the theoretical distributions obtained 
for the normalised scalars for the Gaussian case, we have 
generated 20 fuU-sky CMB Gaussian simulations, in the way 
explained in section 3. A map of the studied normalised 
scalar is then obtained for each simulation and a distribution 
function of the scalar is constructed. Finally, the average 
value and dispersion of the distribution function from the 20 
simulations is computed and plotted versus the theoretical 
distribution. 



4.1 The normalised Hessian matrix scalars 

4.1.1 The normalised eigenvalues 

The normalised eigenvalues are defined in terms of the orig- 
inal eigenvalues through the following expression: 



Al 
A2 



2 ^ 



j_ -I i. 




Therefore analogously to the original eigenvalues, A2 < Ai < 
00 and —00 < A2 < Al. Note that for a power spectrum as 
the one of the CMB, we have a\ « (J2. Using this approxi- 
mation we find for the previous equation that the normalised 
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Normalised scalar Notation Definition Domain 



pdf 



greatest eigenvalue Ai i + (A2,oo) 



p(Ai) 



3V^ 



lowest eigenvalue 



As 



i(A+-A_) (-00, Ai) p(A2): 



^1 1 



^Aie^^ 



1 - er/ I 



Laplacian 



(—00, 00) 



dctcmiiiiant 



shear 



A1A2 {— cxd,cxd) 



y 



(0,00) 



p(d) = 



_4_ 4d 



1 - it! 



d < 

J > 



distortion 


A_ 




(0,00) 




p(A_) = 2A_e"^- 


ellipticity 


e 


eCT2 


(—00, 00) 




p(e) = 4|e| (l + 8e2)"^ 


shape index 


Z 


^arctan(-Jg) 


(-1,1) 






gradient 


9 


9 


(0,00) 






derivative of gradient 






(—00, 00) 






Gaussian curvature 




d 


(-00,00) p{kg) = < 


' -iz 


KG V 3«G 






yi^e'V^e*^ [l-er/(V6i)] dz kg > 



Table 1. List of the normalised scalars. The pdf's are valid for any homogeneous and isotropic Gaussian random field, independently of 
its power spectrum. 



eigenvalues are approximately proportional to the original 
ones and, therefore, have basically the same physical mean- 
ing. Their distribution functions for the Gaussian case are 
given by (see Appendix C): 



p(Ai) = 



3V27r 



4e' 



1 + 



2^ 
Aie 3 



3V27r 



27r. 2- 
-A2e 




The pdf for Ai and A2, together with the corresponding re- 
sults obtained from the simulations, are shown in the top 
panels of Fig. 7. 



If.. 1.2 The normalised Laplacian 

The normalised Laplacian, A+, is defined using the original 
Laplacian: 



0-2 



(61) 



where, by definition —00 < A+ < 00. The normalised Lapla- 
cian of a HIGRF follows a Gaussian distribution with zero 
mean and unit dispersion: 



P(A+) = 



27r 



(62) 



Fig. 7 (bottom left panel) shows that the distribution ob- 
tained for A+ from the Gaussian CMB simulations perfectly 
follows its theoretical pdf. 



J^.l.S The normalised determinant 

The normalised determinant, d, is defined analogously to 
the original one, using the normalised scalars Ai and A2 as 
follows: 



d = A1A2 



(63) 



where d £ (— oo,(X)). Using the approximation ai « a2 
valid for CMB maps, the normalised determinant becomes 
proportional to the original determinant. 

The probability density function of the normalised de- 
terminant for the Gaussian case is (Sec Appendix C). 



p(d) = 



J_ Ad 



1 - er/ 



d < 
cl> 



(64) 



The bottom right panel of Fig. 7 shows that the agreement 
between the theoretical pdf and that obtained from simu- 
lated CMB maps is very good. 
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Figure 7. Theoretical distributions and results obtained from CMB Gaussian simulations for the normalised eigenvalues (top left and 
right), the normalised Laplacian (bottom left) and the normalised determinant (bottom right). The normalised Lapla<;ian is Gaussian 
distributed whereas the other three normalised scalars follow different non-Gaussian distribution. 



4.2 The normalised distortion scalars 

4-2.1 The normalised shear 

The normalised shear, y, is defined as a constant factor times 
the original shear, 



y = 



y 



2ai 



(65) 



so inherited by the original shear, < y < oo. The distribu- 
tion function of this normalised quantity for the Gaussian 
case is simpler than for y : 



p(y) = 4e 



-Ay 



(66) 



4-2.2 The normalised distortion 

The normalised distortion, A_ , is defined by: 

r A- 



(67) 



and analogously to the original distortion, < A- < cxa. The 
probability density function of the normalised distortion for 
a HIGRF is 



p{\-) = 2A-e" 



(68) 



This theoretical distribution versus the one obtained from 
simulations are given in the top right panel of Fig.8. 



The comparison between this probability density function 
and the one obtained from 20 Gaussian CMB simulations, is 
shown in the top left panel of Fig. 8. Note that the agreement 
between the theoretical curve and the results obtained from 
simulations is very good. 



4-2.3 The normalised ellipticity 

The normalised ellipticity, e, is proportional to e : 

(T2 



(69) 
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Figure 8. Theoretical distributions compared to the results from CMB Gaussian simulations for the normalised shear (top left), the 
normalised distortion (top right), the normalised ellipticity (bottom left) and the normalised shape index (bottom right). 



where, by construction —oo < e < oo. The probability den- 
sity function of the normalised ellipticity of a HIGRF is 



p(e) =4|e| (l + 8e'')" 



(70) 



Note that this distribution function has very long tails and 
that < > is not defined. This problem was inherited from 
the original ellipticity pdf, and motivates the introduction 
of the shape index, which is a bounded quantity. 

The theoretical pdf and the results obtained from CMB 
simulations are compared in the bottom left panel of Fig. 8. 



4-2.4 The normalised shape index 

The normalised shape index, Z, is defined in terms of the 
normalised ellipticity in an analogous way to the definition 
of the original shape index, (see equation (28)): 



Z = — arctan 

Tt 



(71) 



Therefore —1 < I < 1. The distribution function of the 
normalised shape index, for the Gaussian case, is deduced 



from the one of the normaUsed ellipticity through simple 

transformations: 



P(0 = ^1 



{¥)\ 



1 + 2 cot' 



(72) 



The bottom right panel of Fig. 8 shows the good agreement 
between the theoretical prediction and the results obtained 
from simulations for the normalised shape index. 



4.3 The normalised gradient related scalars 

4-3.1 The normalised square of the gradient modulus 

The normalised square of the gradient modulus, g, is defined 
as the original gradient times a constant factor: 



(73) 



therefore < g < oo and the resulting distribution function 
of this normalised quantity for the Gaussian case, is simply 



p{g) = e" 



(74) 
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Figure 5. Gradient related scalars associated to tlie temperature 
map of Fig. 2 are shown: the squared modulus of the gradient 
(top), and the gradient derivative (bottom). 




-..Mn. ^^^^^^^^^^^^H -'.1.^-'. 



Figure 6. The maps of curvature scalars corresponding to the 
Gaussian CMB simulation given in Fig. 2 normalised to unity 
dispersion, (so we work with the adimensional map AT/o-g.), are 
shown: Gaussian curvature (top) and extrinsic curvature (bot- 
tom) 



The top panel of Fig. 9 shows the theoretical distribution 
function compared to the results obtained from simulations 
for g, showing excellent agreement. 



4-3.2 The normalised derivative of the square of gradient 
modulus 

The normalised derivative of the square of the gradient mod- 
ulus, Dg, is proportional to the original scalar: 



b - 



\/8cr,.(7p 



(75) 



so inherited by the original scalar, — oo < Dg < oo. The 
probability density function of this normalised scalar for the 
Gaussian case, is trivially obtained from the pdf of the prim- 
itive scalar: 



P{Dg) = 



dy 
y ' 



(76) 



The bottom panel of Fig. 9 shows the theoretical pdf and the 
one obtained from simulations for Dg. Again the agreement 
is very good. 



4.4 The normalised curvature scalars 

4.4-^ The normalised Gaussian curvature 

The normalised Gaussian curvature, ka, is defined in terms 
of the normalised determinant and the normalised gradient, 
as follows: 



d 



KG 



(77) 



where by construction —oo < kg < oo. The pdf of the 
normalised Gaussian curvature of a HIGRF is given by: 



p(kg) = 




[l — erf (A/fe)] dz 



The comparison between the theoretical and simulated re- 
sults are shown in the top panel of Fig. 10. 



4.4-2 The normalised extrinsic curvature 

The normalised extrinsic curvature, kex, is defined analo- 
gously to the original scalar, using equation (35): 



2 VT" 



A-, 



Dg 



(1 + 5) 



(79) 
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Figure 9. Comparison between the theoretical pdf and the one 
obtained from Gaussian CMB siirmlations for the normalised 
square of the gradient modulus (top) and the normalised deriva- 
tive of the gradient. 



Therefore — oo < kex < oo. As for the original Kex, it is very 
complicated to obtain an analytical expression for the pdf 
of this normalised scalar. However we have included in the 
bottom panel of Fig 10, the pdf of the normalised extrinsic 
curvature obtained from 20 Gaussian CMB simulations. 

4.5 Correlations between normalised scalars 

By looking at the definitions of the different scalars, it be- 
comes apparent that some of the scalars arc related to each 
other. In order to know how much independent information 
they contain, it is interesting to study the correlations be- 
tween the different scalars. Using simulations, we have cal- 
culated these correlations for a HIGRF which are given in 
table 2. Wc have included all the normalised scalars except 
for the ellipticity, due to the fact that the dispersion of this 
quantity is not defined and, therefore, the usual correlation 
coefficient can not be calculated. However, the ellipticity and 



8 - 



6 - 



2 - 




-1 -0.5 0.5 




-2-10 1 



Figure 10. Comparison between the theoretical pdf and the 
one obtained from 20 Gaussian CMB simulations for the nor- 
malised Gaussian curvature (top), and distribution function ob- 
tained from simulations for the normalised extrinsic curvature 
(bottom) . 



the shape index are very closely related and thus it is ex- 
pected that the correlations between the ellipticity and the 
rest of the scalars will be similar to those obtained for the 

shape index. 

It is interesting to note that the gradient of the modulus, 
g, is the only scalar which is uncorrelated with all the con- 
sidered scalars. On the contrary, the eigenvalues, Ai and A2, 
have some degree of correlation (or anticorrelation) with all 
the scalars (except for g). Within the Hessian scalars, wc see 
that only the normalised Laplacian, A+, and determinant, d, 
are uncorrelated between them. A+ is also uncorrelated with 
two of the normalised distortion scalars (the shear, y, and 
the distorsion, A_) as well as with Gaussian curvature kg. 
Conversely, d has some correlation with these three scalars 
but is uncorrelated with the normalised shape index, I, the 
derivative of the gradient, Dg and extrinsic curvature, kex- 
With regard to the distortion scalars, y and A- are strongly 
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correlated, which is expected since y is, except by a constant 
factor, the square of A-. However, these two scalars are un- 
correlated with the shape index l. It is also interesting to 
note that kg is uncorrelated with hex- 



5 TREATMENT OF INSTRUMENTAL NOISE 
AND MASKS 

Real data contain not only the cosmological signal but also 
contaminating emissions and instrumental noise. The Galac- 
tic region, where the Galactic foregrounds dominate, is usu- 
ally masked from the data. It is also a common procedure to 
mask the emission coming from extragalactic point sources. 
These masked pixels are then discarded from the analysis. 

In addition, the instrumental noise produces disconti- 
nuities from pixel to pixel in the map, since this noise is 
assumed to be, in general, uncorrelated (white noise). The 
discontinuities introduce problems in the calculation of the 
derivatives. In order to deal with this issue, we will smooth 
the CMB signal plus noise with a Gaussian beam of FWHM 
equal to 2.4 times the pixel size (where we assume that this 
is the size of the Gaussian beam used to filter the original 
CMB signal). 

Assuming that the white noise is homogeneous and 
isotropic, its power spectrum is given by: 



C7 = 



47rcr„ 

^ pixtc 



(80) 



where a„ is the noise dispersion and Npixtot is the total 
number of pixels. So the resulting power spectrum of the 
noisy simulations can be expressed as follows: 



(81) 



where <Tg is the Gaussian beam dispersion and C| is the 
CMB power spectrum. 

However, in most experiments, the instrumental noise is 
usually an anisotropic Gaussian random field, charactarised 
by a different dispersion at each pixel an{x) (with x the uni- 
tary vector on the sphere in the direction of observation). In 
order to deal with anisotropic noise, we have built a "pixel- 
dependent power spectrum" H{{x) in the following way: 



He{x) 



CI 



pixtot 



(82) 



He{x) would be the power spectrum of a map that contain 
the filtered CMB plus isotropic noise with dispersion cr„(x), 
and filtered again with a Gaussian beam of dispcrson ag. In 
practice, due to this second smoothing, the dispersion of the 
noise in a given pixel would depend not only ou the uoisc 
level on that position but also of its neighbours. However, 
provided that the dispersion of the noise varies smoothly, 
which is usually the case, the previous equation is a good 
approximation to the "power spectrum in each pixel" . 

Since, by construction, the normalised scalars are inde- 
pendent of the power spectrum of the field, we can obtain 
these quantities for each pixel, taking into account equation 
(82). Wc just need to construct the moment fields on the 
sphere ao{x),ai{x) and a2{x) using equation (36), where the 
Ci's are now given by H({x). The normalised scalars are then 
calculated by introducing these pixel-dependent moments 



on their corresponding definitions. The probability distri- 
bution function obtained from the map of the normalised 
scalar constructed in this way will follow the corresponding 
theoretical distribution given in the previous section. This 
is one of the advantages of working with the normalised 
scalars. Note that if we construct the ordinary scalars for a 
map containing anisotropic noise, the resulting probability 
distribution will not longer follow the theoretical pdf for a 
Gaussian field. 

In addition to noise wc also have a masked region, which 
introduces discontinuities in the boundary of the mask. 
Moreover, we have no information about the field inside 
these masked regions. As before, in order to deal with this 
problem, the first step is to smooth the masked map (with 
the pixels of the mask set to zero) with a Gaussian beam of 
FWHM equal to 2.4 times the pixel size. As already men- 
tioned, this solves the problem of the discontinuities of the 
noise and also reduces the mask boundary problem. How- 
ever a more sophisticated procedure is necessary to eliminate 
the effect of the mask, since if we calculate the normalised 
scalars from this smoothed masked map, those pixels close 
to the boundary would be strongly contaminated by the 
spurious signal of the mask. Nevertheless, the values of the 
normalised scalars in those pixels far enough from the mask 
would be correct. Therefore we need to generate an extended 
mask that eliminates from the analysis not only the original 
masked pixels but also those pixels in the neighborhood of 
the mask. The particular shape of this extended mask would 
depend on the original mask but also on the particular nor- 
malised scalar that we are considering. 

We can obtain the extended mask for each normalised 
scalars using simulations. In particular, we construct and 
compare one exact and one approximated map for the con- 
sidered quantity for each simulated map. The exact map le 
is constructed, simply, by calculating the normalised scalar 
from the smoothed noisy map; since the derivatives are ob- 
tained without masking any region, they will be properly 
calculated. The approximated map la is obtained by calcu- 
lating the corresponding normalised scalar from the masked 
smoothed noisy map. In this case the derivatives of those 
pixels close to the boundary of the masked region would be 
heavily contaminated by the mask. For each simulation, at 
each pixel outside the original mask, we calculate then the 
error quantity e defined as: 



(83) 



where ae is the dispersion of the exact map for the corre- 
sponding simulation. The average of this quantity over a 
large number of simulations is then obtained and we have 
a map of e at each pixel outside the original mask. The ex- 
tended mask will be formed, in addition to the pixels of the 
original mask, by those pixels with e greater than a fixed 
value e, , that is, we keep for the analysis only those pix- 
els where the value of the considered normalised scalar is 
reasonably close to the correct value and therefore are not 
appreciably contaminated by the presence of the mask. Con- 
sidering only the pixels outside the extended mask, the dis- 
tribution of the normalised scalars should follow the theo- 
retical pdf's found in previous sections. 

In order to test the method using realistic simulations, 
we have studied the normalised scalars using masked CMB 
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0.40 
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0.64 


1.00 


0.91 


0.23 


-0.40 


A+ 


0.91 


0.91 


1.00 
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-0.23 


0.23 





1.00 


-0.58 


y 


0.40 


-0.40 





-0.58 


1.00 


A_ 


0.42 


-0.42 





-0.55 


0.96 


Z 


-0.83 


-0.83 


-0.92 
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Da 


-0.52 


-0.52 


-0.58 








KG 


-0.19 


0.19 





0.83 


-0.48 


Kex 


0.89 


0.89 


0.99 









Table 2. Correlations between the normalised scalars. 



Scalar 


Useful pixels (%) 


Ai 


70.6 


A2 


70.5 


A-H 


70.1 


d 


72.0 


y 


67.8 


X 


()!).2 


I 


68.1 


9 


70.7 


Dg 


71.8 


Hg 


70.7 


Kex 


70.1 



Table 3. The percentage of useful pixels kept after applying the 
extended masks obtained with e* =0.1 for the normalised scalars. 



maps containing anisotropic noise. We have used the KpO 
mask, defined in Bennett et al. 2003b, which covers approx- 
imately a 24 per cent of the sky. The anisotropic noise hcis 
been simulated with a level equal to that expected for the 
noise-weighted average of the Q, V and W WMAP chan- 
nels (Komatsu et al. 2003) after two years of observation. 
The signal-to-noise ratio per pixel of the simulations ranges 
from 1.6 to 5.9 (at resolution Ngide = 256). To obtain the 
extended mask we have used 100 simulations and find that 
a value of e, — 0.1 provides good results. We have then 
obtained from the simulations the pdf for each of the nor- 
malised scalars (but using only those pixels outside the ex- 
tended mask) and compare it with the theoretical pdf. Ta- 
ble 3 gives the percentage of useful pixels that are kept af- 
ter applying the corresponding extended mask for the nor- 
malised scalars. For all the normalised scalars we see that 
approximately a 70 per cent of the pixels can be used for 
the analysis, loosing only around a 6 per cent of the sky 
with respect to the KpO meisk. The extended mask for the 
normalised ellipticity can not be obtained using the error 
quantity e, since the dispersion of this scalar is not defined. 
Alternative error quantities could be used to obtain its ex- 
tended mask. However, for simplicity, we have used the mask 



A- 


I 


g 


Da 

y 


Kq 




0.42 


-0.83 





-0.52 


-0.19 


0.89 


-0.42 


-0.83 





-0.52 


0.19 


0.89 





-0.92 





-0.58 





0.99 


-0.55 











0.83 





0.96 











-0.48 





1.00 











-0.46 








1.00 
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1.00 
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-0.46 











1.00 








-0.91 





-0.57 





1.00 



obtained for the normalised shape index, since these two 
quantities arc closely related. 

The pdf's obtained from 100 siunilations (using only 
those pixels outside the extended mask) for the normalised 
Laplacian, distortion and shape index are compared to their 
corresponding theoretical pdf's showing an excellent agree- 
ment. Similar results are obtained for the rest of normalised 
scalars (including the ellipticity), using their correspond- 
ing extended masks (which are not shown for the sake of 
brevity). 



6 CONCLUSIONS 

In this paper we have introduced several scalar quantities 
on the sphere and we have shown how to calculate them for 

a CMB map directly from the aim harmonic coefficients. 
In particular, we have considered the square modulus of 
the gradient, the Laplacian, the determinant of the nega- 
tive Hessian matrix, the distorsion, the shear, the ellipticity, 
the shape index, the eigenvalues of the negative Hessian ma- 
trix, the Gaussian curvature, the extrinsic curvature and the 
gradient derivative. Assuming a homogeneous and isotropic 
Gaussian field, we have derived analytical or semi-analytical 
expressions for the pdf's of each of the scalars. For conve- 
nience, we define normalised scalars which have more simple 
distribution functions than ordinary scalars. Therefore, for 
these normalised scalars we have tested the theoretical re- 
sults with simulations, showing an excellent agreement. 

We propose these scalars to study the Gaussian char- 
acter of the CMB. In particular, we aim to test their power 
using non-Gaussian CMB simulations based in the Edge- 
worth expansion (Martmez-Gonzalez et al. 2002). Moreover, 
we plan to apply this method in future works to both physi- 
cally motivated non-Gaussian models and the WMAP data. 
Finally, we would like to point out that the study of the 
scalars can be adapted, for instance, to deal with regions or 
extrema above (or below) different thresholds. This will also 
be studied in a future work. 



ACKNOWLEDGMENTS 

The authors thank Patricio Vielva for useful comments. 
CM thanks the Spanish Ministerio de Ciencia y Tec- 
nologi'a (MCyT) for a predoctoral FPI fellowship. RBB 



© 0000 RAS, MNRAS 000, 000-000 



16 Monteserin et al. 



0.4 



0.3 



0.1 








^-202 




- 



12 3 



e 




-1 -0.5 0.5 



Figure 11. The distribution obtained from 100 simulations com- 
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APPENDIX A: DERIVATIVES ON THE 
SPHERE 

In this appendix we derive the harmonic coefficients of the 
first and second derivatives of the temperature field as a 
function of the aem- With these expressions and using the 
HEALPix package it is straightforward to calculate the maps 
of the field derivatives as well as the scalars. 

It is a common procedure to write the temperature field 
on the sphere as a series of harmonic functions: 



(Al) 



The spherical harmonic functions constitute an orthonormal 
complete set. Also, they are eigenvectors of the momentum 
operators and L^: 

LzYem = mYim , (A2) 
L^Yem = e{£ + l)Yem. (A3) 
In addition, the L+ and L_ operators that are defined: 
L±= La:± iLy , (A4) 
act on the spherical harmonic functions in the following way: 

L±Yem = ^/{e±m + l) (^Tm)y^m±i ■ (A5) 

In order to obtain the derivatives of the field, we will write 
them as a function of the previous operators: 

^ = i [e-'^L+ - e'HJ\ , (A6) 

^ = zL.. (A7) 

On the other hand, we will also make use of the following 
common spherical harmonic relations: 



-'•^ sine Ye 



m+l — 



{i-m){i-m + l) 
(2£ + l) (2£ + 3) 

i£ + m){£ + m+l) 
{2i - 1) {2£ + 1) 



Ye+lrn — 



Y, 



l-lm , 



(A8) 



e*** sin 9 Yi„ 



[l + m){l + m + l) 
(21 + 1) [21 + 3) 

{I -m){£- rn + 1) 
{2£ - 1) (2£ + 1) 



Ye+i 



Ye- 



(A9) 



and sin be^ and ctm respectively, as a function of the 
harmonic coefficients acrn of the original field* 



hem = H+ {£, m) ae+i,m + H- {£, m) a^-i,^ , 
where H± are functions of £, m given by 



H+{£,m) = -i£ + 2) 



H- {£, m) = {£- 1) 



(£ + m + l)(^-m + l) 
{2£ + 1) {2£ + 3) 



{£ + m)(£- m) 



(All) 
(A12) 

(A13) 
(A14) 



{2e + 1) (2£ - 1) ' 

Note that these functions are even on m and that H -{£,£) 
is zero for all the £ values. 

Analogously, we can obtain the harmonic coefficients 
dtm, eem, and fem corresponding to the fields sin^ ^^gi" 
and sin 61^: 



-m aen 



(A15) 



eev 



WW 



,m) 



aem + 



+ 
+ 



i£-lf J 
{^)H+{£,m)H+{£+l,m) 



£- 



H- (£, m) H- {£ - 1, m) 



ae+2,m + 
ae-2,m ,(A16) 



im[H+ {£,m)ae+i,m + H- (£,m) ai^i^m] ■ (A17) 



Making use of the previous results and taking into ac- 
count the different expressions presented along this work, 
it is straightforward to obtain the maps of the considered 
scalars. 



APPENDIX B: SOME USEFUL RESULTS 
REGARDING SPHERICAL HARMONIC 
SERIES 

In this appendix we present some useful results for spher- 
ical harmonic series, which are necessary to calculate the 

covariances of the variables defined in equations (5-9). 

In particular, we have made use of the following ex- 
pressions which involve the harmonic functions and their 
derivatives with respect to 9: 



E 



dYem ^* _ „ 



(Bl) 



cos 9 Yem 



{£-m + l){£ + m+l) 
{2£+l){2£ + 3) 



Ye+im + 



{£ + rn) {£ - m) 
{2£ + 1) {2£ - 1) 



Ye- 



(AlO) 



Applying the derivative operators given in (A6) and 
(A7) to the temperature field (Al), and taking into account 
the previous spherical harmonic relations, it is straightfor- 
ward to obtain the harmonic coefficients of the fields S- 



89 89 



9^^. 21 + i i{i + 1) 



47r 



* Note that similar expressions to obtain the first derivatives of 
the temperature field have been independently obtained by Erik- 
sen et al. (2004b). We also want to remark that in Schmalzing 
& Gorski 1998 a similar method to obtain the derivatives on the 
sphere was described. 
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. 



Another interesting set of series involve the presence of 
powers of m. In particular, we need: 



Since S, T and R are independent variables, the numerator 

(B3) 

and the denominator are also independent. 

The numerator is the square root of the addition of the 
square of two Gaussian uncorrelated variables with the same 
dispersion as- Therefore the numerator, z, should follow a 
Rayleigh distribution: 



2 2£+l 

em. 



4-7T 



8 [^(^+1)]'- 2^(^+1) 



cos 261 



,(B4) , . 1 - 



2^+lZ(Z + l) . 2, 

— : — sm 

47r 2 



and applying to the previous expression, we get: 



E 



2dYem^, 21 + ll{l + 1) . „ . 
m „„ Yfrn. = — : — - sm 6 cos t 



(B5) 



(B6) 



It can also be be shown that other series that appear in 
the calculation of the covariances of the variables (5-9) and 
that involve odd powers of m are zero, such as: 



E 

m=-l 

Ei dY, 
m Ye 

m= — 

I 

E 



Im 







m,= -t 

where i is an odd integer. 



(B7) 
(B8) 

(B9) 
(BIO) 



APPENDIX C: NORMALISED SCALARS 
DISTRIBUTION FUNCTIONS 

In this appendix wo give some guidelines on how to deduce 
the pdf of ordinary scalars for a HIGRF. We also show how 
to obtain the probability density function of the normalised 
scalars from that of the original scalars. As an illustration, 
we focus on the normalised ellipticity and determinant, but 
similar calculations are needed to obtain the pdf 's of the 
rest of the scalars. 

As explained in section 2, ordinary scalars can be con- 
structed in terms of the Gaussian variables {p, g, r, s, t} given 
in equations (5) to (9). In particular, the ellipticity can be 
rewritten as: 



2(r + s) 



(CI) 



Let us define the uncorrelated Gaussian variables 7? = r + s, 
S = r — s and T = 2t. It is straightforward to show that 
aji = ffi and ag — = — af. In terms of these new 
variables, the ellipticity is given by 

-I- T2 



2R 



(C2) 



(C3) 



where z > 0. The denominator is Gaussian distributed with 
zero mean and 2aR dispersion. Performing a change of vari- 
able we can obtain the pdf of the inverse of the denominator 

X = 



1 

2R- 



(C4) 



where — oo < a; < oo. Taking into account that x and z are 
independent, we have: 



p{x,z) =p{x)p{z) = 



s e 



(C5) 



Since we need to calculate the pdf of e = xz, we perform a 
change of variables (a:, z) (e, z) and then, integrating over 
z (note that z > 0), we obtain the distribution function of 
the ellipticity given in equation (49). 

The normalised ellipticity is proportional to the ordi- 
nary ellipticity, e = ce, and, therefore, the pdf of e can be 
trivially obtained from p{e) by a simple change of variables: 



p(e) = lp(e=^) 



(C6) 



The preceding expresion leads to equation (70). The same 
applies for all the normalised scalars which are proportional 
to the ordinary ones. 

However, some normalised scalars are related to the 
ordinary ones in a more complex way. To construct their 
pdf's it is convenient, whenever possible, to rewrite them 
in terms of A+ and A_, since these two normalised scalars 
arc independent. This can be proved using the property 
PAi(*) = PA2(— a;), already mentioned in section 3, where 
p\^ denotes the pdf of \i. As an example, we show how to 
calculate the probability density function of the normalised 
determinant, d, which can be expressed as: 



d = A1A2 



Al 



(C7) 



Using the probability distribution function of A+ and A-, 
(equations 62) and 68) and by performing different changes 
of variables, we can obtain the distribution function of the 
normalised determinant. First we calculate the pdf of the 
u = X\ and v = XL, which are given by: 



p{u) 
p{v) 



1 



2n \fu 



e 2 



where < u < 00 and < < 00. Since A+ and A- 
independent, u and v are also independent, therefore 



p(u,v) =p(u)p{v) 



1 



27r a/m 



1 

e e ■ 



(C8) 

(C9) 
_ are 

(CIO) 



We perform then a change of variables (u, v) — > (w, d), where 
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d = J — J. Finally, integrating over u we obtain the pdf of 
the normalised determinant d given in equation (64). 

The normalised eigenvalues can also be rewritten in 
terms of A+ and A_, through the following expression: 

(t)-(i4)(fc)' 

Therefore, throught a straightforward change of variables, 
the pdf of the normalised eigenvalues can be obtained from 
those of A+ and A- . 

The pdf 's of the remaining normalised scalars can be 
constructed in an analogous way, since they can be trivially 
written cither in terms of A+ and A- or of the other nor- 
malised scalars. 
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